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Abstract 

We investigate a model containing two species of one-dimensional fermions 
interacting via a gauge field determined by the positions of all particles of the 
opposite species. The model can be solved exactly via a simple unitary trans- 
formation. Nevertheless, correlation functions exhibit nontrivial interaction- 
dependent exponents. A similar model defined on a lattice is introduced and 
solved. Various generalizations, e.g. to the case of internal symmetries of the 
fermions, are discussed. The present treatment also clarifies certain aspects 

of Luttinger's original solution of the "Luttinger model" . 
71.10.-w, 71.27.+a 
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Exactly solvable models Jl|-§] have played an important role in the current understand- 
ing of one-dimensional interacting many-particle systems. These together with the idea 



gence of the "Luttinger liquid" 0] as a unifying concept. Nevertheless, the technicalities of 
these exact solutions (bosonization, Bethe ansatz) are often rather complex. In the present 
paper we wish to introduce a class of interacting models which can be diagonalized by a 
simple (pseudo-)unitary transformation, yet exhibit nontrivial Luttinger-liquid behavior. 
The models can be defined both in the continuum and on a lattice, and can have rather 
arbitrary single-particle bandstructure. Only the interactions are constrained to be of a par- 
ticular "gauge form". The long-distance asymptotics of correlation functions can then be 
determined either by direct calculation or by a mapping on an effective Luttinger liquid de- 
scription. Our investigation was inspired by Luttinger's original treatment of the "Luttinger 
model" , and we will comment on this connection further below. 

We start by considering the simplest model in our class, a one-dimensional fermion model 
with two species of particles, designated by a pseudospin index a = ±, having coordinates 
x a i and momenta p a i = —id Xai . The Hamiltonian of our model then is 



where we have introduced a "covariant momentum" IT CT j = p^i + aA^^Xai), i.e. in this model, 
particles interact via a gauge potential, given for a particle at x by A a (x) = J2j V(x — X- a j)- 
The potential V is an even function and vanishes at infinity. On a ring of length L, we will 
assume a potential V that is periodic in L. Clearly, the Hamiltonian is not time-reversal 
invariant, but it is invariant under simultaneous time reversal and charge (a) conjugation. 

The model can now be straightforwardly diagonalized by a (in general pseudo-)unitary 
transformation: noting that 



of dominant low energy bosonic excitations of Fermi systems H gave rise to the emer 




at 



(1) 



e iS({x +i },{x-j}) 



iS{{x +i },{x-j}) 



p +i - d x+i S({x +i }, {x-j}) 



(2) 
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one can chose S so as to eliminate the interaction in eq.(|l|) by 

S({x +i },{x- j })=^2E(x +i -x- j ) , (3) 
where E is the indefinite integral of the interaction potential: 

E(x) = f X dx' V(x') . (4) 



o 

The transformed Hamiltonian then takes the form 

H = j s He~ iS = Wpli ■ (5) 

The eigenfunctions of H clearly are Slater determinants of plane wave states ({A;-^}, {k + j}), 
characterized by the sets of wavenumbers {k- t i} and {k + j} for the — and + particles re- 
spectively Consequently, the eigenfunctions and eigenvalues of the original Hamiltonian are 
obtained straightforwardly: 

He- iS \{k^},{k +J }) = \Y.kl % e- lS \{k.,},{k + ,}) ■ (6) 

cri 

At first sight it thus appears that the spectrum of the interacting Hamiltonian is independent 
of the interaction. Conformal field theory, or equivalently Luttinger liquid theory, then 
would imply that the asymptotic form of correlation functions (which is directly determined 
by the eigenvalue spectrum) is also interaction-independent. This conclusion is however 
incorrect: periodic boundary condition have to be treated carefully. In fact, keeping all 
other coordinates fixed, one easily finds S{x-^ = L) — S(x-j = 0) = — N + S, where L is the 
length over which periodic boundary conditions are applied, N + is the total number of + 
particles, and the phase shift 5 is given by 

1 dxV(x) (7) 

An analogous result, with iV + and 5 replaced by iV_ and —5, holds for the phase shift of the 
+ particles. Consequently, the quantization condition on the wavenumbers is given by 

Lk±j =f N T S = 27rn± t j , (8) 
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where the n±^ are integer quantum numbers analogous to those used in the noninteracting 
case. Clearly, particles of one given "spin" orientation give rise to an effective Aharonov- 
Bohm flux acting on the other species, the value of the flux depending on the number of 
particles present. It should now also be clear why we refer to the transformation eq.(0) as 
pseudoumtaiy: unless 6 is "accidentally" an integer multiple of 2n, the plane wave states of 
the interacting and noninteracting problems obey different boundary conditions and there- 
fore define different Hilbert spaces. 

The ground state energy Eq can be found in any sector with N± particles, as follows. In 
order to minimize E we must choose n± t i = n° ± i ± [[^riV T ]] irii , where n® ± i are the quantum 
numbers in the absence of the interaction, and we denote x = [[x]]i nt + [[x]) rem for any x 
where is the closest integer to x. Thus —1/2 < [[x]] refn < 1/2 for any x. The change 

in energy due to the interaction, 5E = ^r{iV+[[^]]r em + ^-[[^~]]rem}) * s n °t extensive, 
but on the scale expected from a magnetic field applied to the ring. 

An effective Luttinger liquid description in terms of a bosonic field theory for the low- 
energy properties can be obtained from the low-energy excited states. |J To be precise, 
we start from a ground state with N± = 2n + 1 and assume that N± Q 5 is an integer 
multiple of 2tt. We now add n±R (u±l) particles at the right (left) Fermi points of the ± 
particles. Introducing particle number and current quantum numbers N± = n±n + n±i and 
J± = n±R — n±L the second order variation of the ground state energy is 

£ {2) = ^(2^o + 1) [(tt 2 + 5 2 )(Nl + N 2 ) + vr 2 (J 2 + J 2 ) + 2n5( J + iV_ - J_iV + )] .(9) 

Up to quantum fluctuations, N± and J± are related to bosonic fields and their conjugate 
momentum density via N± = —(L/ix)d x (f)± and J± = ZTI-i-. The effective Hamiltonian 
including the low-energy quantum fluctuations then takes the form 



H = - f dx { [l + (5/tt) 2 ] [{d x <P + f + (d x <f> 



7T 2 



ni + n 2 
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+2S\n.d x <f> + -u + d x <i>.]} , (io) 
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where n = 4no/L is the particle density. Introducing new variables 

0± = 0± , tl± = U±T- 2 d x ct>T (11) 

the Hamiltonian takes an apparently non-interacting form (eq.(^) with 5 = 0). However, 
the expression of single-fermion operators |J is changed and therefore the asymptotic decay 
law of the single particle Green function is obtained as 

G R ±(x) = (V« ± (x)V>k(0)> w e^ ±N ^ x x- l ~ a , (12) 

with a = S 2 /(2tt 2 ). Thus, for any non-vanishing 5 the decay is faster than 1/x, leading, 
amongst other things, to the well-known power-law singularity of the momentum distri- 
bution function at kp. The correctness of eq.(|T2"D can be checked independently using the 
eigenfunctions of eq.(|6|): one obtains a Toeplitz determinant of the form previously consid- 
ered by Luttinger ||, and which has the same asymptotic power law as obtained by the 
bosonization approach. 

Similarly, correlations of two-particle operators decay as x~ v , with interaction dependent 



exponent r/. Specifically: 

VWl± ^i = 2 (13) 

4^2* V2 = 1 + (1tS/tt) 2 (14) 

ipR±ipL± r/ 3 = 2 + 2(5/vr) 2 (15) 

?Pr±?Pl t => m = 1 + (1 ± 5/n) 2 (16) 



The most slowly decaying correlations identify the dominant incipient instabilities. In the 
spin language, for positive 5 then spiral spin-density wave correlations and opposite-spin 
Cooper pairing correlations with one fixed spin orientation (jj. and || are not degenerate) 
are favored, whereas for negative 5 correlations with reversed spin orientations dominate. 
Adding a density-density interaction between the two spin orientations, the degeneracy be- 
tween pairing and spin-density wave correlations is lifted. The density correlations, eq.(|T^), 
are not affected by the interactions because they are diagonal element of the density matrix 
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which themselves are unchanged by the unitary transformation, eq.(0). We notice that the 
exponent for pairing correlations with equal pseudospin, eq. (p3|) , is just twice the exponent 
of the single-particle Green function, i.e. there are no singular vertex corrections in this 
particular two-particle correlation function. Finally, from eq.(|8]) it is clear that the value 
of 5 is relevant only modulo 2tt. Consequently, the results (|T^) to flTfip are valid only for 
\8\ < 7T. Outside this interval 5 has to be taken modulo 2tt. We note that the scaling rela- 



tions between the different exponents in eqs. ( \T@i ) to ( \Tq ) are different from those of standard 



fermionic Luttinger liquids because of the presence of time-reversal breaking terms in the 
Hamiltonian. 

We can now comment on Luttinger's original solution of his model. In first-quantized 
form his Hamiltonian only contains first derivatives: 

#Lu tt = E ^ = E^-E p-j + 2 E v(*-3 ~ x +i) • ( 17 ) 

ai i j ij 

We first remark that this Hamiltonian is a conserved quantity as far as the Hamiltonian ([!]) 
is conserved, and shares all (non degenerate) eigenfunctions. It is unbounded from below 
though, unlike H in Eq (]]]). Hence the issue of finding its groundstate is replete with diffi- 
culties familiar from relativistic field theories. The second-quantized version of the model 
-f^Lutt can be solved consistently and exactly by filling the Dirac sea and using bosoniza- 
tion. || This leads, amongst other things, to an asymptotic decay exponent of the single 



particle Green function a = 1/yl — (S/tv) 2 — 1. In a first-quantized framework, a consis- 
tent but different solution can be obtained if one is willing to consider quasi-groundstates 
where single-particle states below a certain very negative cutoff energy -E cu toff are left empty 
(evidently, the model does not have a conventional groundstate). This "rapidity cutoff" in 
fact is frequently used in the Bethe ansatz solution of field theoretical models Such 
a state becomes natural if one is interested in finding the groundstate of H in eq.(|l|), and 
examines the eigenvalue of i?Lutt> a commuting operator, in this state. The transformations 
used for eq.(^) also can be used here, and the solution found as earlier and lead to a shift 
in its eigenvalue due to interactions SE Lutt = jjiV+K^]],™ + AL[[^^]] rem }, a number of 



the 0(1) as one would expect from a current carrying state. The correlation function can 
be found using Luttinger's original paper and lead to the same asymptotic decay exponent 
a of the Green function as in eq.flT^p. This result was in fact obtained in Luttinger's paper, 
i.e. Luttinger's result in fact applies to the first-quantized solution of the model described 
here. |5j The same result for correlation exponents can also be obtained by considering vari- 
ations of the energy with particle number, similar to what we described above. We note 
that the Mattis-Lieb and Luttinger results for a, though different in general, agree to the 
lowest nontrivial order in S. The differences at higher order clearly have to be attributed to 
the different cutoff procedures used in the two calculations. 

We now turn to similar models defined on a one-dimensional lattice. Specifically, we will 
consider the Hamiltonian 

# = -EE 

a m=l 

where a is a periodic function, ot m+ L = a m , and the number operator n m = cl n c m . This 
corresponds to a lattice where the up electrons feel a "gauge potential" due to the down 
electrons, and vice versa. Here periodic boundary conditions are implied, i.e. cl+i,o = c\ a . 
We notice that for the particular case where only «o is nonzero, this lattice model only 
involves two-particle interactions, contrary to our continuum model where some three-body 
interaction is unavoidable. 

The interaction term V can be either zero, or one of two non trivial functions which 
retain the exact solvability. We will consider the XXZ model and the Hubbard model, i.e. 

Vxxz = V^2n ja n j+ha and VH ub = U ^ %T n iJ. ( 19 ) 
j> j 

and will show below that these are exactly solvable. The XXZ model corresponds to two 
copies of the usual model, wherein the two species of particles (spin up and down) only talk 
to each other via the phase factors. The Hubbard model corresponds to the usual two body 
interaction. 

We now perform a unitary transformation induced by U = exp(iS) where S = 
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exp io 



E 

z 



a. 



-in. 



-TO+1<T 



+ h.c. 



+ V 



(18) 



J2i<i, m <L A,m^zT n »4- ^ ^ S eas y t° see that (3ij = —fijj , i.e. an odd function is appro- 
priate, and we will assume it to be so. Thus we find 



Cmcr * 6 C mr7 6 C ma exp 



-ia ^2 Pmj.ni,- 
i 



(20) 



The transformed Hamiltonian takes the form 

~L-1 / 



eX P ia ^{^m,l ~ An+1,2 + OW-jK-" C mer 



.m=l 



- exp ^wr^{/? Lii - + ax-i}ni,-<rj ^ L(T c Xa + h.c. 



(21) 



We now use the freedom in defining (3 to cancel the interior terms in the phase factor by 
choosing 



Pm+i,i ~ Pm,l = a m-i for 1 < m < L - 1 and 1 < I < L. 
The hop across the L «-> 1 bond has a total phase 

I 



(22) 



(23) 



It is in fact not necessary to solve explicitly for (3, although it is easy enough to do so for 
simple choices of a. We can add eqs. (p2[) for 1 < m < L — 1 and further add to it to 
find 

L 



- Pl,l + ®L-l = a «" 
n=l 
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(24) 



since the sum is independent of I. This gives \ a = N_ a a5. The number operator N a — ► iV^ in 
any sector, and hence we see that the problem collapses to one with lattice fermions having 
twisted boundary conditions. If V = 0, we can follow the logic used for the continuum 
model to determine asymptotics of correlation functions. It turns out that up to the trivial 
replacement vf = 7iN/2L — > 2sm(nN/2L) one obtains the same expression for as in 
the continuum limit, and consequently the same low-energy effective Hamiltonian, eq.(|lO|). 



and the same expressions for correlation exponents (eqs.flEj) to (16)) apply. 
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In the presence of a nonzero extra interaction V, previous work |T(J can be used where 
the Bethe Ansatz has been adapted to the case of a "spin twist" , which is precisely the case 
needed here. We write the solution immediately: for the XXZ model 

£=-£cosA£ (25) 

n,cr 

and 

LK = 2ni: + aSN-. + £ 0(A£ - C) (26) 

m 

with the usual phase shift and the usual integers I n . For the Hubbard model we have a pair 
of equations. With E = — cos k n with 1 < n < N, and N — + we find 

Lk n = 2nl n + N l 5 + 2J2 tan" 1 [4(A,- - sin k n )/U] 

j 

2£tan~ 1 [4(A j -smk n )/U] = 2nJj - 5N + 2 £ tan" 1 [2(A j - A k )/U]. (27) 

n k 

In these models, one has non Fermi liquid behavior even in the absence of a r , and adding this 
changes the exponents, and indeed even the symmetries of the model, e.g. for the Hubbard 
model we have less than SU(2) invariance. The detailed behavior of the lattice models and 
the resulting exponents will be reported elsewhere. 

There is a number of further possible generalizations of the present model: first, the 
unitary transformation (fj) will in fact put Hamiltonians containing arbitrary powers of the 
covariant momenta H a i into diagonal form. Provided the highest nonvanishing power is 
even, these models have a well-defined groundstate. One could thus study models with 
complicated bandstructures, involving e.g. more than two Fermi points. Similarly, in the 
lattice model certain forms of hopping terms beyond nearest neighbors can be included. 

Another generalization is obtained by giving additional internal degrees of freedom to 
the a = + and a = — particles. For example, assuming that they both occur in m different 
"flavors" one obtains a model with an internal SU(m) x SU(m) symmetry. By a calculation 
analogous to that leading to the exponents in eqs.(|13D to (|T6|) one finds r/i = 2, i] 2 = 2 =p 
25/n+mS 2 /n 2 = 2+2m5 2 /n 2 , rj^ = 2±25/ir+m5 2 /TT 2 . As expected from symmetry, these 
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exponents are independent of the flavor indices appearing in the corresponding operators. 
One can further solve the case where the number of flavors for the + and — particles is 
different. 

In conclusion we have presented a class of lattice and continuum fermion models which 
are exactly solvable by a pseudo-unitary transformation, leading to nontrivial and non- 
Fermi-liquid behavior, with exponents depending upon the interaction. The models, unlike 
those solvable by bosonization, do not have an unbounded spectrum, and eliminate the 
problem of the negative energy Dirac sea and consequent Schwinger terms, and help us to 
focus on the physics of the interactions in one dimension in a bounded, and even a finite 
dimensional Hilbert space (for the lattice models). The method used embeds the original 
problem considered by Luttinger in a family of commuting Hamiltonians which contain 
both bounded as well as unbounded operators. By focusing on the problem of finding the 
groundstate of the bounded operators one comes up with eigenfunctions which are of the 
type considered by Luttinger, enabling us to make a connection between the methods used 
by him (Toeplitz determinants and the Szego formula for asymptotics) with more recent 
conformal/Luttinger liquid methods. The relatively simple form of the exact wavefunctions 
may also make it possible to understand in detail physical properties in the non-asymptotic 
(intermediate and high energy) regime where many questions still remain open, even in 
otherwise well-understood one-dimensional models. 

One of us (H.J.S.) wishes to acknowledge the warm hospitality of the Indian Institute of 
Science, Bangalore, where some of this work was done. This research was supported in part 
by the National Science Foundation under Grant No. PHY94-07194. 
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